Non-entropic theory of rubber elasticity: 
flexible chains with weak excluded- volume interactions 



A.D. Drozdov* 
Department of Chemical Engineering 
Ben-Gurion University of the Negev 
P.O. Box 653 
Beer-Sheva 84105, Israel 



Abstract 

Strain energy density is calculated for a network of flexible chains with weak excluded- volume 
interactions (whose energy is small compared with thermal energy) . Constitutive equations are 
developed for an incompressible network of chains with segment interactions at finite deforma- 
tions. These relations are applied to the study of uniaxial and equi-biaxial tension (compression), 
where the stress-strain diagrams are analyzed numerically. It is demonstrated that intra-chain 
interactions (i) cause an increase in the Young's modulus of the network and (ii) induce the 
growth of stresses (compared to an appropriate network of Gaussian chains), which becomes 
substantial at relatively large elongation ratios. The effect of excluded-volume interactions on 
the elastic response strongly depends on the deformation mode, in particular, it is more pro- 
nounced at equi-biaxial tension than at uniaxial elongation. 

Key-words: Flexible chain, Excluded-volume interaction, Polymer network. Finite deformation. 
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1 Introduction 

This study is concerned with the elastic response of polymer networks at finite strains. According 
to the classical theory of rubber elasticity [Tj, an arbitrary chain in a network is treated as Gaus- 
sian, which allows a simple formula to be derived for the strain energy density of a network, and 
stress-strain relations to be developed in the analytical form. Two shortcomings of the concept 
of Gaussian chains are traditionally emphasized: (i) this model does not account for long-range 
interactions between segments (Gaussian chains can intersect themselves), and (ii) the end-to-end 
distance of a Gaussian chain exceeds its contour length with a non-zero probability. To avoid these 
disadvantages, it seems enticing to replace Gaussian chains in a network by fiexible chains with 
excluded-volume interactions (this model does not permit self-intersections) or by semi-flexible 
chains (this approach guarantees that the end-to-end distance of a chain is always less than its 
contour length). Although the necessity to go beyond Gaussian chains has been realized for a long 
time (the seminal paper by Flory |2] appeared more than half a century ago), our knowledge of 
the mechanical behavior of networks of non-Gaussian chains remains rather limited, due to some 
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difficulties in their mathematical treatment. Serious progress in the analysis of statistics of sin- 
gle polymer chains and membranes with excluded-volume interactions was reached by using the 
renormalization group technique |31ISE1IS1- These methods, however, have been employed for the 
analysis of the distribution functions only and have not yet been applied to determine the strain 
energy of a network, despite the importance of the latter problem for applications, see |H1 and 
the references therein, as well as recent publications [51 [TUl [TTl [T^ . 

This study is motivated by a problem which, at first glance, appears to be quite simple. Consider 
an incompressible permanent network of Gaussian chains with a given set of parameters (segment 
length bo, contour length L, number of chains per unit volume M) at a fixed absolute temperature 
T. Suppose that Gaussian chains in the network are replaced by flexible chains (with the same 
parameters bo, L and M) with excluded-volume interactions (whose strength vq will be defined 
later). The question is how the presence of intra-chain interactions affects the Young's modulus 
E of the network? In a more general context, this question may be reformulated as what is the 
influence of segment interactions on the stress-strain relations for the network? 

In order to shed some light on this issue, it is necessary to define unambiguously what flexible 
chains with excluded-volume interactions mean. Two approaches are conventionally employed to 
describe configurations of a polymer chain. According to the first, a chain is thought of as a random 
walk with a small step length bo and a large number of steps A'^ (their product L = Nbo is assumed to 
be finite when — > oo). Excluded-volume interactions between segments are treated as a constraint 
that rules out trajectories that intersect themselves (self-avoiding random walks), which implies that 
the "strength" of segment interactions vq has no physical meaning. A disadvantage of this concept 
is that the distribution function for end-to-end vectors Q of self-avoiding chains is unknown. Some 
approximations for this function are available at Q ^ oo [7:, but they are insufficient for the 
determination of elastic moduli. Nevertheless, the effect of intra-chain interactions on the elastic 
moduli may be assessed qualitatively: as the number of configurations for a self-avoiding walk is 
lower than that for a walk without constraints, and the free energy of a chain is proportional to its 
entropy (which, in turn, is proportional to the logarithm of the number of available configurations), 
it is plausible to assume that excluded-volume interactions between segments reduce the stiffness. 

According to the other approach, a chain is treated as a curve with length L in a three- 
dimensional space. Any configuration of the chain is described by the equation r = r(s), where r 
stands for the radius vector, and s G [0, L]. This configuration is characterized by a weight (energy), 
which is determined by some Hamiltonian H{r). For a Gaussian chain, the Hamiltonian reads 

while for a chain with excluded-volume interactions, this functional is given by ^3] 

H{r) = Ho{r) + ^r) (2) 

with 

m = ^ J^"^ J\{r{s)-r{s'))dsds'. (3) 

Here ks is Boltzmann's constant, 6{r) denotes the Dirac delta-function, and the pre-factor vq 
characterizes strength of segment interactions. 

Unlike the concept of random walks, an assessment of the influence of segment interactions 
on the Young's modulus E becomes non-trivial in this case. On the one hand, the presence of 
the second term on the right-hand side on Eq. (2) reduces the number of configurations with a 
noticeable weight, which results in a decrease in the free energy. On the other hand, this term 
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increases the energy H of any available configuration, which causes the growth of the average 
energy of a chain, and, as a consequence, an increase in the strain energy density of a network. 

A rigorous treatment of the interplay between these two factors is the objective of this paper. 
Due to some technical difficulties in the evaluation of path integrals for the Hamiltonians (2) and 
(3), we confine ourselves to the analysis of weak excluded-volume interactions, whose energy $(r) 
is small compared with thermal energy ksT. 

The exposition is organized as follows. Section 2 deals with the free energy of a flexible chain 
and the strain energy density of a permanent network of polymer chains calculated within the 
entropic and non-entropic concepts. The free energy of a flexible chain with a small, but arbitrary 
functional $(r) is found in Section 3 in terms of an appropriate correlation function. Section 4 
has a merely technical character. We derive some explicit expressions for correlation functions, 
which are employed in Section 5 to determine the free energy of a chain with excluded-volume 
interactions. Constitutive equations for a network of chains with weak segment interactions are 
developed in Section 6, where they are applied to the analysis of uniaxial and equi-biaxial tension 
of an incompressible medium. Some concluding remarks are formulated in Section 7. 



2 The concept of non-entropic elasticity 

We begin with a brief exposition of the classical theory of rubber elasticity demonstrate its 
shortcomings, and introduce some refinement of the conventional approach. The concept of entropic 
elasticity is grounded on the treatment of a polymer chain as a random walk in a three-dimensional 
space. For definiteness, we suppose that the walk begins at the origin and has a fixed length 6o of 
each step. The main hypothesis of this theory is that the distribution function p(Q) of end-to-end 
vectors Q entirely describes configurations of a chain. The free energy of a chain ^'(Q) is connected 
with the distribution function p(Q) by the Boltzmann equation 

p(Q) = exp(-|M). (4, 

In the nonlinear elasticity theory, two states of a medium are distinguished: (i) the reference 
(initial) state occupied before application of external loads, and (ii) the actual (deformed) state 
that is acquired after deformation. As a polymer chain is entirely characterized by the relative 
positions of its end-points, two vectors are introduced: the end-to-end vector in the reference state 
Q, and that in the actual state Q'. These quantities obey the equality 

Q' = F • Q, (5) 

where F is a deformation gradient, and the dot stands for inner product. It follows from Eqs. (4) 
and (5) that the increment of free energy 

AM/(F,Q) = vI/(Q')-^(Q) 

driven by deformation of the chain reads 

A^'(F, Q) = kBT\\np{Ct) - lnp(F • Q) . 

The strain energy per chain VF(F) is determined by averaging the increment of free energy over 
the initial distribution of end-to-end vectors, 



W{¥) = ksT J [inp(Q) - lnp(F • Q)]p(Q)fiQ, (6) 
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where the integration is performed over the entire space. Given a strain energy W, stress-strain 
relations for a chain are determined by conventional formulas, see, e.g., |14| . Equation (6) is 
noticeably simplified when the distribution function p(Q) is isotropic: 

p{Ci)=p*{Q), Q = \Q\- 

Bearing in mind that |F • Q| = (Q • C • Q)2, where 

C = F 



is the right Cauchy-Green deformation tensor, and T stands for transpose, and introducing a 
spherical coordinate frame {Q, </>, 9}, we find that 



r2n 

-i2. 



W{C) = kBT p^{Q)Q^dQ d(t> [Inp^Q) - lnp,(^(Q • C • Q)2 jj smOdO. (7) 
In particular, for a Gaussian chain with the radial distribution function 

'"<« = (2l?)'-K-§)- (8) 

where b = ^/boL is the mean square end-to-end distance, Eq. (7) implies the classical formula 

W{C) = hBT{liiC)-3), (9) 

where Im stands for the mth principal invariant of a tensor. 

According to the other way of modeling a polymer chain each configuration is associated 
with a curve r(s) in a three-dimensional space. For a chain that begins at the origin and finishes 
at a point Q, the radius vector r(s) satisfies the boundary conditions 

r(0) = 0, r(L) = Q. (10) 

The Green function (propagator) of a chain whose energy is described by a Hamiltonian H{r) reads 

nr{L)=Q / H(r(s))\ 

G{Q)= exp U-ZZ Vp[r(s)], (U) 



r(0)=0 V ^B^ 

where the path integral with the measure P[r] is calculated over all curves r(s) that obey Eq. (10). 
As the functional integral is determined up to an arbitrary multiplier ^S] , the additional restriction 
is imposed on the function G(Q), 

'"G(Q)dQ = l, (12) 



which ensures that the Green function G(Q) coincides with the distribution function of end-to-end 
vectors p(Q). 

Within the entropic elasticity theory, the strain energy W oi a flexible chain with a Hamiltonian 
H is determined by using the same technique as for a chain treated as a random walk: given H{r), 
the Green function G(Q) is calculated from Eq. (11) and is normalized with the help of Eq. (12) 
to obtain the distribution function p(Q). Afterwards, the strain energy T4^(F) is determined by Eq. 
(6). 

Two shortcomings of this approach should be mentioned: (i) for a chain with a Hamiltonian H, 
a correct measure of the free energy is the average Hamiltonian, while the use of Eqs. (6), (11) and 
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(12) appears to be unnecessary and questionable, and (ii) our previous analysis of a flexible chain 
grafted on a rigid surface demonstrates that the entropic elasticity theory leads to conclusions that 
contradict physical intuition pT)| . 

In this study, we associate the free energy ^(Q) of a chain with an end-to-end vector Q with 
the average Hamiltonian of this chain 



§(Q) 



{H)q 
1 



G(Q) 7r(0)=0 



r(L)=Q 



H{r{s)) exp 



Hjrjs)) 



V[r{s)]. 



(13) 



Given ^(Q), we determine the normalized end-to-end distribution function p(Q) from the equation 
similar to Eq. (4), 



P{Q) = exp 



^(Q) 



exp 



(14) 



and calculate the strain energy VF(F) as the average (over the distribution function) increment of 
free energy. 



W{F) 



§(F.Q)-§(Q) p(Q)dQ 



(15) 



Simple algebra reveals that for a Gaussian chain with Hamiltonian (1), our approach coincides with 
the conventional one and results in the strain energy density (9). However, for a flexible chain with 
excluded-volume interactions, Eqs. (13) to (15) differ from appropriate relations developed within 
the entropic elasticity theory. 

An important remark regarding Eqs. (13) to (15) is that an additive constant in the expression 
for the free energy ^'(Q) does not affect the strain energy density M^(F) (as it is expected). This 
constant is excluded from the formula for the distribution function p{Q) by the normalization 
condition (14), whereas it disappears in Eq. (15) because only the increment of free energy is 
substantial for the determination of VF(F). 



3 Free energy of a chain with weak segment interactions 



Our aim now is to calculate the free energy of a flexible chain with weak intra-chain interactions, 
Eqs. (2) and (3), in the first approximation with respect to the ratio VQ/ksT. Inserting expressions 
(2) and (3) into Eq. (11), expanding the exponent into the Taylor series, and disregarding terms 
beyond the first order of smallness, we obtain 



G(Q) 



r{L)=Q 
r(0)=0 



1 



keT 



exp 



kBT 



V[ris)] 



1 



Q 



ksT 



Go(Q) 



(16) 



where Go(Q) is the Green function for a Gaussian chain, and the superscript index zero stands for 
averaging with the help of the Gaussian Hamiltonian. It follows from Eqs. (2), (3) and (13) that 
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G(Q) 7r(0)=0 



r(L)=Q 



Fo(r(s)) + $(r(s)) exp 



1 



kBT 



{Ho{r{s)) + ^r{s)) 



V[ris)]. 



Expanding the exponent into the Taylor series, using Eq. (16), and neglecting terms beyond the 
first order of smallness, we find that 



Q 



1 



ksT / Go(Q) 7r(o)=o 



r(L)=Q 



{Ho{r{s)) + $(r(.)) - J-Ho{r{s))^T{s)) 



kBT 



5 



X exp 



1 



Hoir{s)) 



V[ris)] 



(17) 



Our purpose now is to determine all terms on the right-hand side of Eq. (17) separately. To find the 
average of the Gaussian Hamiltonian {Ho)q, we insert expression (1) into Eq. (11) and introduce 
an explicit dependence of the Green function Gq on segment length bo 



fr{L)=Q r 3 fL ^j^j. 2 



rr(L)=Q 
r(0): 

Differentiation of Eq. (18) with respect to 60 implies that 
dG, 



V[r{s)]. 



(18) 



3 r'^fdr n2 ■ 



1 fr{L)=Q 



ksTho 7r(o)=o 



-ffo(r(s))exp 



2^0 J a ^ds 
Hoiris)) 



P[r(.)] 



V[r{s)] 



1 



ksTbo 

It follows from this relation that 



(i/o)^Go(6o,Q). 



/u \0 _ ksTbo dGp 

^""'^"^-G^M)^^''^^- 

Substitution of expression (8) with b^ = boL into Eq. (19) yields 



(19) 



_ ^ksT ^Q^ ^ 



(20) 



According to the definition of ($)q and Eq. (1), we have 

rr(L)=Q r 3 fL^^j^^. 

$(r(s))exp -- 

r(0)= 



/•r(L)=Q r 3 /-i/dr \2 ■ 



where an explicit dependence of (^)q on bo is omitted for brevity. Differentiation of this equality 
with respect to 60 results in 



^(W2,Go((«,Q)) 



/:::r(Mr(S('))^-)*« 



r(0) 

X exp 



s)) 
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1 /•r(L)=Q 



V[ris)] 



ksTbo 7r(o)=o 



i7o(r(s))$(r(s))exp 



Hoirjs)) 
ksT 



V[r{s)] 



(21) 
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It follows from Eq. (21) that 



Go{bo,Q)dbo 



(($)^Go(6o,Q)) 



ksTbo 



~dbr, 



ksTbo dGp 
QGo(6o,Q) dbo 



(&o,Q). 



Combination of Eqs. (19) and (22) implies that 



Inserting expressions (20) and (23) into Eq. (17), we arrive at the formula 



(22) 



(23) 



(24) 



According to Eq. (24), to find the free energy of a flexible chain with weak segment interactions, 
we need to determine the energy of interactions averaged with respect to the Gaussian Hamiltonian 
($)q only. To calculate this quantity, a perturbative technique is applied. 



4 Correlation functions 

Our aim now is to derive an explicit expression for the average of an exponential function of r(s) 
when the averaging is performed with respect to the Hamiltonian Hq. We begin with the perturbed 
Green function 



cm) 



r(t)=Q 
r{0)=0 



exp 



— k^r~ Jo ^^^^ ' ~ ^ 



V[ris)] 



= (exp(^^p(s)-r'(s)d5))^Go(Q), 



where p(s) is a smooth vector function determined on [0, L], 

r°(5) = q| 

is the "classical" path for a Gaussian chain, and 

r'(s) = r(s)-r°(s). 
Substitution of Eqs. (26) and (27) into Eqs. (11) and (25) results in 

2fe2 



G^(Q) = exp(-^)GPQ), Go(Q) = exp(- 



G{Q), 



where 



exp 



GP(Q) = 



V[r'{s)], 



V[r'{s)], 



(25) 



(26) 



(27) 



(28) 



(29) 
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and the path integrals are calculated over all curves r'(s) that satisfy the boundary conditions 

r'(0) = r'(L) = 0. (30) 

Setting 

r'(.s) = ro(s) + R(s), (31) 
where the functions ro(s) and R(s) obey the zero boundary conditions 

ro(0) = ro(L) = 0, R(0) = R(L) = 0. (32) 
we transform the expression in the square brackets in the first equality in Eq. (29) as follows: 

3 r^/dr' 



_3 

'2^0 JO 

rL 



+ 



/ fp(s) • ro(s) + p(s) • R(s) 



S 



c/s. 



(33) 



Integration by parts with the use of Eq. (32) yields 

dro , . dR, 



ds'^ 



{s)ds. 



It follows from this equality and Eq. (33) that 



A = , 

26o Jo 

r-L 

+ 



ds + 



3 d^R 



6o ds^ 



(s)+p(s) •ro(s)ds 



/ p(s) • R(s)ds. 
io 



Assuming the function R(s) to satisfy the differential equation 

3 d^R, 



bo ds^ 



we find that 



A 



3 f^wdro /dR. 

Jo y'^^^'^J + i^^^ 



(s) + p(s) = 0, 



26o 



ds + p(s) • R(s)ds. 



The solution of Eq. (34) with boundary conditions (32) reads 

bo '•^ 



where 



Kis) = fJ^ Dis,s')pis')ds', 
Dis,s') = s'{l-j) is>s'), Dis,s')=s{l-j) is<s'). 



(34) 



(35) 



(36) 



Integrating by parts the second term in the square brackets in Eq. (35) and using Eqs. (32), (34) 
and (36), we obtain 

d^R, 



ds''' 



'-{s)ds 



bo 



p(s) • R(s)(is 



= {jf D{s, s'Ms) ■ pis')dsds'. 
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Substitution of this expression and Eq. (36) into Eq. (35) results in 



dri 



L rL 



(38) 



Bearing in mind that the last term in expression (38) can be taken away from the path integral, 
we conclude from Eqs. (28), (29) and (38) that 

"6n /-^ /-^ 



GP(Q) = Go(Q)exp 



^ / / D{s,s')^p{s)■p{s')dsds' 
o Jo Jo 



Comparison of this equality with Eq. (25) implies that 

rL . rL rL 



r f 1 ^0 f f 

(exp[y^ p(s) • r'(s)dsJ)Q = exp — D{s, s')p{s) ■ p{s')dsds' 



(39) 



Equation (39) serves as the basic tool for the analysis of correlations between the radius vectors 
r(s) at various points s G [0, L] of a chain. Setting 

p(s) = qS{s - ti) - qS{s - t2), 

where q is an arbitrary vector, and ti,t2 G (0, L) are arbitrary points, we find from Eq. (39) that 



(exp[q • (r'ih) - r'(i2))])^ = exp(^A(ti, ^2)) , 



6 

where 

A(ii,i2) = D{h,h) - 2D(tut2)+D{t2M)- 
In particular, when q = — zk, where i = \/— T, and k is a real vector, Eq. (40) reads 



(40) 
(41) 



>Q = exp 



6 



(exp[-ik- {r'{ti) -Y'{t2 
Replacing r' by r in accord with Eqs. (26) and (27), we obtain 

(exp[-zk • (r(i2) - r(ti))])^ = (exp[-zk • (r'{h) - r'{t2 



A(ii,i2) . 



exp 



(-^A(ii,t2)-zk.Q 



^°Qexp(-^k.Q^) 

(42) 



6 -^"-"-^ - L 

Our aim now is to apply Eq. (42) in order to calculate the average energy of segment interactions 
($)q for the functional $(r) given by Eq. (3). 

5 The average energy of segment interactions 

We introduce the Fourier transform of the delta-function by the formula 

<5(r) = ^ / exp(-^k-r)a!k, (43) 

combine Eqs. (3) and (43), and find that in a spherical coordinate frame {A;, ^, ^}, whose 63 vector 
is directed along the end-to-end vector Q, the functional <I> reads 



$(r) = n k^dk d(f) r sinOde ds exp[-ik • fr(s) - r(s' 

2L^(27r)'^ Jo Jo Jo Jo Jo ^ 



vo 





/ k'^dk / d(j) / sinOde / ds / exp -ik - (r(s)-r(s')) 
Jo Jo Jo Jo Jo ^ ' 



ds' 
ds'. 
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It follows from this relation and Eq. (42) that 



L2(27r)3 7o 
x/%xp(- 



roo rliv PIT pL 

/ k'^dk / # / sin^d^ / ds 

Jo Jo Jo Jo 



■^A(s, s') - ikQ cos e'—-^)ds'. 



Equations (37) and (41) imply that for any s > s', 

Using this equality, performing integration over 0, and setting t = s — s', we obtain 



^^0 



(27rL)2 
(27rL)2 Jo 



roo rn rL rs p bok'^ / \ T 

/ k'^dk / sin 0^6* / ds exp (r — — ) — ikQ cos 6— 

Jo Jo Jo JO L 6 V L / L 



k^dk I sinede 
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Jo 



dT 



T) exp 



bok^ 



{r-'-)-^kQ 



cos 9 



L 



dr, 



where we changed the order of integration over s and r and integrated over s explicitly. Introducing 
the notation x = cos 6 and t = t/L, we find that 



exp 



{bkf 



t{l -t)- ikQxt (1 - t)dt. 



Bearing in mind that 



we arrive at the formula 



1 , ,^ 2sin(A;Qt) 

exp(—ikQxt)dx = — — , 

-1 kQt 



^ ^ 
27r2 Jo 



Mo{kQ,k)k^dk, 



(44) 



where 



and 



Mo{a, k) = / m{at) exp -^^t{l - t) (1 - t)dt 
Jo L 6 -I 



m{x) = 



sma; 



X 



Differentiation of Eq. (44) with respect to bo implies that 



6o- 



dbo 



Vo 
"27r2 



where 



Mi{a,k) 



m 

6 



2 rl 



m{at) exp 



/•oo 

/ Mi{kQ,k)k^dk, 
Jo 

{bkf 



(45) 



6 



t{i - 1) t{i - tydt. 



Substitution of expressions (44) and (45) into Eq. (24) yields 



m) = ^-^{^-l) + ^,j^ M{kQ, k)k'dk, 



(46) 



where 



M(a, k) = m{at) exp 
Jo 



^^til-t)](l+^^t{l-t)\{l-t)dt. 
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Changing the order of integration in Eq. (46) , we obtain 



where 



J{t,Q) 



exp 



62 
{bkf 



til-t) 



27r2g Jo t 

{bkf 



J{t, Q)dt, 



(47) 



1 + 



t{l - t) sm{kQt)kdk. 



Bearing in mind that the function under the integral is even, we present this equahty in the form 



J{t,Q) = \ I exp 



{bkf 
~6~ 



t{l-t) 



1 + i^f(l - t) ) sm{kQt)kdk. 



The integral in Eq. (48) is calculated with the help of the formulas 



/ exp( —)sm{l3k}kdk = \ exp -— , 

J-oo ^ 2, ' V Q: Q: \ Za^ 



-oo 



J^J-p{-^)^Mmkdk = ^\/^^exp(-£-)(3 - 



1 /2^/3 
2 



that are fulfilled for any a > 0. Combination of these relations with Eq. (48) implies that 



J(t,Q) 



3QH ^ / 



3QH 



i). 



262(1 -i) 



Inserting this expression into Eq. (47), we find that 



*(Q) = 

where v = vo/{kBT), and 



SksT 



(91 
V62 



1 + 



4^ 



(5i?i(g) 



RiiQ) 
R2{Q) 



= Jo ^"K-2P(rr^) 



62 



dt 



R2{Q 



exp 



(- 



262(i-t);^i3(i_i)' 

3Q2t ^ dt 



262(1 -t)Vi(l-i)^' 



The function R2{Q) reads 



26 



exp 



_6 [2^ 

Q 



(48) 



(49) 



(50) 



(51) 



where we used the following variables: t = t/{l — t), s = -y/r, and z = Qs\/3/b. The first integral 
in Eq. (50) is presented in the form 



Ri{Q)=Ro- £ l-exp(- 



3QH 



262(1 -i). 



dt 



(52) 



where Ro is independent of Q. According to the remark at the end of Section 2, the additive 
constant Ro does not affect the increment of free energy A^* and the distribution function p, and 
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we do not calculate this quantity. Some concern may arise regarding Rq, because the integral in 
the formula for Rq diverges. This divergence does not affect, however, the free energy. To avoid it, 
one can replace Eq. (43) for the delta-function by its regularization, 



S{r) = J [/(fe)exp(-zk.r)dk, 



where 



U{k) = 1 iO<k< h), U{k) = 0, {k> K) 



with some fc* ^ 1, find the free energy of a flexible chain with the regularized potential of segment 
interactions, disregard the additive constant in the expression for ^'(Q), and, afterwards, take the 
limit at fe* ^ (X). We do not dwell on detailed transformations, because they cause unnecessary 
complications of the analysis without influence on the final result. 
The term dependent on Q in Eq. (52) is transformed as follows: 



/■If / ^QH ^] dt r-r / 3QVx 



Q 
b 



262(1 -i) 

(1 — exp(— s)Js 2(isj 

where we set r = t/(l — t) and s = '^Q'^t / (2lP'). Integration by parts results in 

(l - exp(-s))s-ids = 2j^ y-^ eM-y)dy = 2r(-) = 20F, 
where T{x) is the Euler gamma-function. Combining these relations, we find that 

RiiQ) = Rq-jVQ^. 



T ^dr 



(53) 



Substituting expressions (51) and (53) into Eq. (49) and neglecting additive constants, we arrive 
at the formula 



(54) 



where 



9w 
7r6^ 



is the dimensionless strength of excluded volume interactions. It follows from Eqs. (14) and (54) 
that 

^2 g^ 



P(Q) = Pq exp 



1(91 



2 V 62 5 , 

where the pre-factor pQ is determined from the normalization condition 

1 r /- r 3/g2 Qxi 2 ,^1-1 

""'^^XL ^"Pi-2(^-^-6)J^^^/ • 
For weak excluded- volume interactions (e <C 1), the integral in Eq. (56) is given by 



(55) 



(56) 



exp 



2V62 



Q 



Q'^dQ = 



3\/3 

6^ 
3^3 



exp 



z^dz + 



:V3 
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where z = Q-y/^/h. This relation together with Eq. (56) imphes that 



(57) 



in the first approximation with respect to e. Inserting Eqs. (54) and (55) into Eq. (15) and 
introducing a spherical coordinate frame {Q,^, we obtain 

r2-K 



w = 



exp 



3/Q 
2U2 



b 



Q'dQ 



/^^ (Q • C • Q - Q2) - ^ (v/QTcTQ - q) 



sin^d^. 



(58) 



Denote by the eigenvectors of the the right Cauchy-Green tensor C and by \m appropriate 
eigenvalues. If the 63 vector of the spherical coordinate frame is directed along the eigenvector is, 
the expression Q • C • Q reads 

Q ■ C • Q = Q'^S, S = (Ai cos^ (f) + X2 sin^ (p) sin^ 6* + A3 cos^ 9. 

Substitution of this relation into Eq. (58) results in 

W = Wi- W2, 



(59) 



where 



Wi 



W2 











/ exp 






zUz I 


Jo 






Jo 



2tt 



J ((Ai cos^ (f) + X2 sin^ 0) sin^ 6* + A3 cos^ ^) ^ ^ ^^^^ 



SeksTpob^ 



f 

Jo 



exp 



z-^dz 



j-2-K 

/ dct> 
Jo 



((Ai cos^ ^ + A2 sin^ (f)) sin^ 6' + A3 cos^ 6") ^ - 1 sin ed0, 



and z = Q/b. We calculate the integrals over z, substitute expression (57) for the coefficient po, 
disregards terms beyond the first order of smallness with respect to e, and find that 

sin 6d9, 



Wi = f 1 + Iv -) / ^ # r ff (^1 cos^ + ^2 sin^ (p) sin^ 6 + X3 cos^ o) - 1 

87r ^ 3 V TT' 70 io '-^ ' 



Wo = 



^^^^ y^^^" # [((Ai cos^ + A2 sin^ 0) sin^ ^ + A3 cos^ o)^ - 1 



Ait 



sin Ode. 



(60) 



Calculating the integrals over (f) and 9 in the first equality in Eq. (60), we arrive at the formula 

= ^(1 + 1^1) (Ai + A2 + As - 3). (61) 

According to Eqs. (9) and (61), at £ = (no intra-chain interactions), Wi coincides with the strain 
energy of a Gaussian chain. Setting x = cos 9 in the other equality in Eq. (60) and using the 
evenness of the function under the integral, we find that 



Wo 



eksT [6 

27r V TT 

6 



2tt 



= eksT 



IT 

eksTi - 

V TT 




1 

2 



dx 



d(p J ((Ai cos^ 4> + \2 sin^ ^)(1 - x^) + 

1 . , 1 ■ 

Ai cos^ ^ + A2 sin^ (?!)) (1 — a;^) + X^x^ jdx — 1 



27r 



TT Jo Jo ^ 



Ai cos^ ^ + A2 sin^ — x^) + Xsx^^ ^ dx — 1 



(62) 
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Although it is possible to develop an analytical expression for the strain energy W2 for an arbitrary 
three-dimensional deformation, an appropriate formula is rather cumbersome, and we do not present 
it for the sake of brevity. We confine ourselves to a particular case of axisymmetric deformation 
with 

Ai = A2 = A (63) 

for two reasons: (i) under condition (63) the governing relations remain relatively simple, and (ii) 
deformation processes (63) are typical for experiments on uniaxial and equi-biaxial extension of 
elastomers. Performing integration over (p in Eq. (62), we obtain 



W2 = eksT 

Calculation of the integral over x implies that 

A 



1 

X + {X3- X)x'^)^dx-1 



W2 

W2 
W2 



ekeT 6 
2 

eksT [6 



A3 + 



VA 



In 



A3 + 



3-A 
A 



A3 + VA3 — A 

Vx 



'A- A3 



arcsm - 



ekBT^[^s-l) (A3 = A). 



(A3 > A), 
(A3 < A), 



(64) 



To obtain the strain energy W2 as a function of the principal stretches Am , the parameter A in Eq. 
(64) should be replaced by i(Ai + A2). 



6 Constitutive equations for a network of flexible chains 



To develop stress-strain relations for a permanent network of flexible chains with excluded volume 
interactions, we adopt the following hypotheses: (i) the motion of chains is affine, which means 
that the deformation gradient F coincides with the deformation gradient for macro-deformation, 
and (ii) inter-chain interactions are accounted for by using the incompressibility condition, which 
implies that the strain energy of a network equals the sum of strain energies of individual chains 
|13j . Denote by M the number of chains per unit volume. It follows from Eqs. (59), (61) and (62) 
that the strain energy density (per unit volume of the network) reads 



W{X, 



ksTM r 



-2e 



-r 

vr Jo 







1 , . 1 ■ 

(Ai cos^ (/) + A2 sin^ 0)(1 - x^) + X3X^)^dx - 1 



The principal Cauchy stresses are expressed in terms of the strain energy W as |14| 

dW 



Am 



dXr, 



(65) 



(66) 



where P stands for pressure. Formulas (65) and (66) provide the stress-strain relations for an 
incompressible network of flexible chains with weak excluded-volume interactions. 

Our aim now is to apply these equations in order to evaluate the effect of segment interactions 
on the elastic response of a polymer network at uniaxial tension (compression) and equi-biaxial 
tension. 
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6.1 Unicixial tension 

Uniaxial tension of an incompressible medium is described by the formulas 

xi = k~^Xi, X2 = k~^X2, X3 = kXs, 

where {Xjn} and {xm} are Cartesian coordinates in the reference and actual states, respectively, 
and k denotes elongation ratio. The right Cauchy-Green tensor C is given by 

C = A;~-^(eiei + 6262) + ^^6363, 

where are base vectors of the Cartesian frame in the initial state, and its eigenvalues read 

Ai = A2 = k-\ A3 = A;2. (67) 
Inserting expressions (67) into Eq. (65) and using Eq. (64), we find that 



W = 
W = 



ksTM 
2 

ksTM 



k' + 



0/ ln(^/PWp— T)_ 
7rr+ VA:(P - 1) 



„\ /o/, arcsin \/l — k^ 
3) -e\ -ik + , / ... - 2 



(k > 1), 
{k < 1). (68) 



According to Eq. (66), at uniaxial tension (compression) the longitudinal Cauchy stress S is 
determined as 

^ . dW . dW 8W IdW 
Substitution of expressions (67) into this equality yields 

which implies that the engineering tensile stress a = T,/k is given by 

dW 
a = . 

dk 



(69) 



Combination of Eqs. (68) and (69) implies that 



ksTM 



( 



3Vfc^(A;^ - 1) - (4A;3 - 1) ln(VP + VW^) 



1 + 



{k > 1), 



a 



ksTM 



e 6 



1 

fc2 



s 6 



X 1 



3^/W(T^^^ + (1 - 4/c3) arcsin VT^H^ 



{k < 1). 



(70) 



It is easy to check that lim^^i cr{k) = 0, which means that the reference state is stress-free. 
Differentiation of the first equality in Eq. (70) with respect to k yields 



da 
dk 



3e [6 



8k2 V vr 



(A;3-l)l 



{8k^ + 1) ln(VP + Vk^ - 1) - (lOk^ - l)Jk^{k3 - 1 



(71) 
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The Young's modulus of a network is defined as 

da 



k=l 



Inserting expression (71) into this equaUty and applying the L'Hospital rule to calculate the limit 
at — 1, we find that 

where 

£^0 = SksTM 

stands for the Young's modulus of a network of Gaussian chains. According to Eq. (72), excluded- 
volume interactions induce an increase in the elastic modulus of a network of flexible chains which 
is proportional to the dimensionless strength of segment interactions e. 

To assess the effect of intra-chain interactions on the stress-strain diagrams, we perform numer- 
ical simulation of Eqs. (70) at uniaxial tension and compression. The results of numerical analysis 
are presented in Figures 1 (tension) and 2 (compression), where the reduced tensile stress 



k-k- 



is depicted versus elongation ratio k (the Mooney-Rivlin plots). Here o"o = f^/iksTM) stands for 
the dimensionless engineering tensile stress. In these figures, the functions a^{k) for a network of 
Gaussian chains are presented by horizontal lines, whereas appropriate dependencies for a network 
of flexible chains with segment interactions demonstrate a monotonic increase in tensile and com- 
pressive stresses (it is worth noting that a weak minimum of cr^, at compression observed in Figure 
2 does not reflect non-monotonicity of the dependence cj(/c); the latter function decreases with k in 
the entire domain G (0, 1)). The results plotted in Figures 1 and 2 indicate that excluded- volume 
interactions cause the growth of stiffness of a polymer network, in agreement with Eq. (72) that 
determines the elastic modulus at small strains. 

Our aim now is to assess the influence of segment interactions on the mechanical response of 
an incompressible polymer network at equi-biaxial tension. 

6.2 Equi-bictxial tension 

Equi-biaxial tension of an incompressible material is described by the formulas 

xi = kXi, X2 = kX2, = fc~^X3, 

where k stands for elongation ratio. The right Cauchy-Green tensor reads 

C = A;^(eiei -I- 6262) -I- ^"^6363, 

and its eigenvalues are given by 

Ai = A2 = k^, A3 = k-\ (73) 

As equi-biaxial tests on an incompressible layer are conventionally performed in the tensile mode, 
we confine ourselves to the case A: > 1. Substituting expressions (73) into Eqs. (65) and using Eq. 
(64), we find that 



I e /6"n/ ,2 1 o\ /6"/ 1 k'^ . V/e^ - 1 



k^ 



(74) 
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According to Eqs. (66), the Cauchy tensile stress S is determined by the formula 



S = 



1 dW 



As 



dW kdW 



d\3 2 dk ' 



It follows from this relation that the engineering tensile stress a = T,/k reads 



IdW 
2~dk' 



(75) 



Substitution of expression (74) into Eq. (75) results in 



a = ksTM 




6 r 2 fc3 . ^6 _ 4 




ttU^ k<^-i\yWzri 



arcsm 



It follows from Eq. (76) that limj^^i a{k) = 0, which means that the initial state is stress-free. 

The dependence of the dimensionless tensile stress ctq on elongation ratio k is depicted in Figure 
3 for £ = (a network of Gaussian chains) and £ = 0.5 (a network of flexible chains with excluded- 
volume interactions). This figure shows that segment interactions induce an increase in the tensile 
stress at all elongation ratios k. The difference between tensile stresses monotonically grows with 
A;. At A; = 6.0 (which is in the range of deformations reached in experiments on elastomers), the 
engineering stress in a network of chains with segment interactions exceeds that in a network of 
Gaussian chains by 28%. 

The results of numerical simulation at uniaxial tension of an incompressible medium are also 
presented in Figure 3 for comparison. According to this figure, at relatively large elongation ra- 
tios {k > 3.0), the difference between the tensile stresses (corresponding to these two deformation 
modes) in a Gaussian network disappears, whereas an appropriate difference in a network of flexible 
chains with excluded-volume interactions remains quite pronounced. Although under both defor- 
mation programs, excluded-volume interactions induce an increase in stiffness of the network, the 
effect of segment interactions on the stress-strain relation is stronger at equi-biaxial tension than 
at uniaxial tension. 

7 Concluding remarks 

The concept of non-entropic elasticity is proposed for the analysis of the mechanical response of a 
network of flexible chains with excluded-volume interactions. Unlike the classical theory of entropic 
elasticity, where the free energy of a chain is entirely characterized by the distribution function of 
end-to-end vectors, we associate the free energy of a chain with the average value of its Hamiltonian. 
The average free energy of a chain and the strain energy density of a network are calculated explicitly 
(under the assumption that the strength of excluded-volume interactions is small compared with 
thermal energy) . Constitutive equations are developed for a network of macromolecules with intra- 
chain interactions under an arbitrary three-dimensional deformation. These relations are simplified 
for uniaxial tension (compression) and equi-biaxial tension of an incompressible medium at finite 
strains. An explicit expression is derived for the clastic modulus of a network of chains with weak 
segment interactions. The effect of intra-chain interactions on the stress-strain diagram is analyzed 
numerically. It is demonstrated that excluded-volume interactions result in an increase in the tensile 
stress at both deformation modes under consideration. This growth is substantial (about 30% at 
elongations typical for experiments on rubbers), which means that the account for excluded-volume 
interactions is quite important for applications. 
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List of figures 

Figure 1: The dimensionless reduced stress cr* versus elongation ratio k at uniaxial tension. 

Figure 2: The dimensionless reduced stress cr* versus elongation ratio k at uniaxial compression. 

Figure 3: The dimensionless tensile stress ctq versus elongation ratio k at uniaxial (fihed circles) 
and equi-biaxial (solid lines) deformations. 
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